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Adiabatic rapid passage in quantum 
dots: phonon-assisted decoherence 
and biexciton generation 



(N 

o 

(N 

> 

o 
O 

m 

(N 



K. Gawarecki 12 , S. Luker 2 , D. E. Reiter 2 , T. Kuhn 2 , M. Glassl 3 , V. M. Axt 3 , A. Grodecka-Grad 4 , P. 
Machnikowski 1 * 

1 Institute of Physics, Wroclaw University of Technology, 50-370 Wroclaw, Poland 

2 Institut fiir Festkorperfheorie, Universitat Miinster, Wilhelm-Klemm-Str. 10, 48149 Munster, Germany 

3 Theoretische Physik III, Universitat Bayreuth, 95440 Bayreuth 

4 Niels Bohr Institute, University of Copenhagen, 2100 Copenhagen 0, Denmark 

Received XXXX, revised XXXX, accepted XXXX 
Published online XXXX 



CD 



Key words: quantum dot, adiabatic rapid passage, optical control, correlation expansion, time-convolutionless, path integrals, biexciton 
generation 



Corresponding author: e-mail Pawel.Machnikowski@pwr.wroc.pl 



I 

O 
o 



We study the evolution of a quantum dot controlled by a 
frequency-swept (chirped), linearly polarized laser pulse 
in the presence of carrier-phonon coupling. The final 
occupation of the exciton state is limited both due to 
phonon-induced transitions between the adiabatic spec- 
tral branches and because of phonon-assisted transitions 
to the biexciton state. 



When the biexciton shift is large enough, the quantum 
dot can be modeled as a two-level system, which corre- 
sponds to excitation with circularly polarized light. For 
this case, we compare different methods of simulations: 
(i) a time convolutionless method, (ii) correlation expan- 
sion and (iii) path integrals. We show that results ob- 
tained from these methods agree perfectly at low tem- 
peratures. 
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1 Introduction The interest in the optical control of 
the states of a quantum dot (QD) has been stimulated by 
recent experiments, where the QD state was controlled ei- 
ther by means of Rabi oscillations [ 1 , 2 1 or by using chirped 
laser pulses exploiting the adiabatic rapid passage (ARP) 
OS) . The usage of chirped pulses offers the possibility of a 
quantum control, which is not sensitive to small variations 
of the pulse area as Rabi oscillations are. In both cases, the 
role of phonons has been shown to be crucial, because their 
interaction with the QD leads to a reduction of the con- 
trol efficiency (5j6]|7|. For the ARR it was shown that at 
low temperatures the influence of phonons depends on the 
sign of the chirp. For negative chirps the control is limited 
by phonon interaction while for positive chirps no damp- 
ing is seen for an excitation with circularly polarized light 
l6l . When under linearly polarized excitation a phonon- 
assisted transition to the biexciton is possible thus hinders 
the control of the exciton state also for positive chirps [ 8J. 
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To describe the phonon influence on the ARP, many 
theoretical approaches have been used like correlation 
expansion 07), time-convolutionless approach (SJ, path 
integrals and the Bloch-Redfield-Wangsness theory J9)- 
Though the methods differ in their complexity, they 
all show a phonon-induced damping. In this paper, we 
want to make a quantitative comparison between three of 
those methods: (i) correlation expansion (CE), (ii) time- 
convolutionless method (TCL) and (iii) path integrals (PI). 
We furthermore discuss when the limiting case of a two- 
level model is a good assumption. 

2 Model We consider a QD modeled as a four-level 
system driven by a laser pulse with either circular or linear 
polarization. The Hamiltonian of the excitonic system is 



H x = E Q (\<t+)(<t+\ + k_X<r-|) + (2£o + E B )\B)(B\. 
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|0) denotes the empty dot, \a±) are the two exciton states 
with different circular polarizations, \B) represents the 
biexciton state, Eq is the single-exciton transition energy 
and Eb is the biexciton shift. The states with circular po- 
larization |cr±) are related to the 'linearly polarized' states 
\X),\Y) by \a±) = (\X) ± i\Y))/y/2. A Gaussian laser 
pulse with initial pulse area and initial pulse duration To 
is frequency modulated by a chirp filter with coefficient a. 
The envelope of the chirped pulse is then 

with the duration r = (a 2 /r§ + Tq ) 1 / 2 and the frequency 
chirp rate a = a /(a 2 + Tq) (6). The central frequency 
of the laser pulse luq is tuned to the polaron shifted single 
exciton resonance. 

In the rotating wave approximation, the Hamiltonian 
describing a carrier-light interaction for a + circular polar- 
ization takes the form 

Bgr> = ^ GOXo + l + W-){B\) e-^ 2 /2 + hx . 

Because the ground state is coupled only to |o+), the sys- 
tem can be reduced to a two-level model. For X linear po- 
larization the Hamiltonian reads 

Ht n) = ^ (|0X*l + WW) e^ 2 + h.c. 

Here the coupling is 1 0) <-> \X) <-> \B). The other state 
(\Y)) is uncoupled, therefore only three-levels are needed 
iflOl . Finally, in order to simplify the numerical simula- 
tions, a unitary transformation to the non-uniformly rotat- 
ing frame is made (SJ. 

For the phonon coupling only longitudinal acoustic 
phonons are taken into account. For the typical parame- 
ter range used in experiments EE), the generation of op- 
tical phonon modes is unlikely and can be neglected. The 
carrier-phonon Hamiltonian is H[ nt = S <£> R, where 

S=\<7 + )(a + \ + \a_)(a_\+2\B)(B\ 

and 

R = ^g k b k +h.c. 
fe 

g k are the coupling matrix elements (8} and bk are the 
phonon annihilation operators. The energy of free phonons 
is described by Hp^ = J^k ^hOjfihi where ujk = ck are 
the phonon frequencies. We consider a spherical QD with 
a size of 5 nm. All of the parameters have been taken from 
0. 

We perform the simulations using three different meth- 
ods. Because all of the methods are well described in the 
given references, we just state shortly the main aspect of 
the different approaches: (i) In the TCL method, which 
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Figure 1 (Color online) Final occupation of the single- 
exciton state as a function of initial pulse area at three 
different temperatures: (a) T = 1 K, (b) T = 20 K and 
(c) T = 50 K. Red lines show results for a = ps 2 , green 
ones for a = —20 ps 2 . Dashed lines are data obtained from 
the TCL, solid lines represent CE and points show results 
from PI. 



we implement in the 2nd order in the phonon coupling, 
only the reduced density matrix is represented in the sim- 
ulation and, in addition, no explicit memory structure is 
present 01 Illl2l l8l. (ii) The correlation expansion contains 
higher order density matrices, in particular, phonon related 
variables 1 13 14|. We include terms up to 4th order in the 
phonon coupling, (iii) The PI method directly calculates 
the time evolution operator at discrete times 01511161 . The 
obtained solution is numerically exact and thus provides a 
benchmark for approximate methods like the TCL or the 
CE. For the Rabi oscillations we already showed that the 
CE and PI agree perfectly for sufficiently low temperatures 

03. 

3 Results We start the discussion with the two-level 
model, where the QD is excited by circularly polarized 
light. The simulated evolution of the exciton occupation 
obtained by the three different methods for two differ- 
ent chirps a = —20, ps 2 is shown in Fig. Q] For zero 
chirp we find the well known Rabi oscillations, which are 
damped by phonon interaction. For a = —20 ps 2 the ARP 
leads to an occupation for the exciton state, which is also 
damped by phonon interaction. For increasing temperature 
the damping becomes stronger. 

For both values of a at T = IK, the results show 
an excellent agreement of the different simulation meth- 
ods (Fig. QI a ))- At a moderate temperature of T = 20 K 
the accordance between the TCL and the CE is also rea- 
sonable as shown in Fig. Q2b). Note, that the mismatch 
between the TCL and the CE is larger for Rabi oscilla- 
tions (q = ps 2 ) than for the ARP (a = -20 ps 2 ). 
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Figure 2 (Color online) Absolute value of the difference 
between the final occupation of the single-exciton state cal- 
culated using TCL and CE method as a function of the 
initial pulse area and the chirp a at: (a) T = 1 K, (b) 
T = 20 K and (c) T = 50 K. 



At T = 50 K the mismatch is getting higher (Fig. QI C ))- 
For the Rabi oscillations a renormalization of the oscilla- 
tion period due to the carrier-phonon interaction occurs. 
This is not well described in the TCL, such that the oscil- 
lation periods obtained by the two methods do not agree 
anymore leading to the observed difference. For the ARP 
at a = —20 ps 2 , these renormalizations are not such im- 
portant and the results agree reasonably well even for this 
temperature. These approximate methods can be tested by 
comparing the results to the data obtained by PI. We find 
that results obtained from the CE and the PI agree perfectly 
for all temperatures considered in this paper. 

In order to investigate quantitatively the mismatch be- 
tween the TCL and the CE methods, we calculated the ab- 
solute value of their difference in a wide range of param- 
eters. Fig. |2] presents these results as a map. As it can be 
seen in Fig. |2ja), at T = IK, the absolute value of the 
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Figure 3 (Color online) Final occupation of the single- 
exciton state as a function of the initial pulse area for 
different binding energies Eb at T = 1 K at three differ- 
ent chirp values: (a) a = —20 ps 2 , (b) a = ps 2 and (c) 
a = 20 ps 2 . 



difference is below 3 %. Furthermore, in this case, the mis- 
match is higher for negative values of chirp. The reason is, 
that at low temperatures, for negative values of the chirp 
the correction to the final occupation from carrier-phonon 
interaction is much higher than for positive one JSJ. As 
shown in Fig. [Hb), the discrepancy at T = 20 K can be 
up to 15% for the small chirp a. Fig. [Tfc) shows, that at 
high temperatures (T = 50 K), near the Rabi oscillation 
region, the mismatch between methods can reach 30%. As 
already seen in Fig.QJc), the large difference between TCL 
and CE is mostly related to the difference in oscillation pe- 
riods for the Rabi oscillations. For chirped excitation this 
renormalization does not play a big role, such that here the 
difference between the two methods is rather small. 

In the next step, we consider a linearly polarized laser 
pulse. We investigate the influence of the biexciton shift 
Eg on the decoherence of the single-exciton state. We in- 
crease the value of Eb successively from Eb = — 1 meV to 
Eb = —4 meV and compare these results with the ones ob- 
tained from the two-level system (for circularly polarized 
light), which should agree for very high biexciton shifts 
Eb 0. The final occupation as a function of the pulse area 
for three different chirp values is presented in Fig. [3] For 
a = —20 ps 2 in Fig.[3{a) the phonon coupling to the exci- 
ton state dominates the damping of the ARP and phonon- 
induced transitions to the biexciton state play a minor role 
even for the very small biexciton shift Eb = — 1 me V. In 
consequence, the discrepancy of the single-exciton occu- 
pation between the two- and four-level system is small. 
Fig. |3jb) shows that in the case of Rabi oscillations at 
o = ps 2 for small biexciton shifts an additional damping 
related to phonon-assisted biexciton transitions appears. In 
particular, for high pulse areas a strong deviation from the 
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two-level model is observed. For a = 20 ps 2 shown in 
Fig. [3]c) no dephasing due to phonons takes place in the 
case of the two-level system. Here, for the small biexciton 
shifts Eb = — 1 meV and Eb = —2 meV, the phonon- 
assisted biexciton generation leads to a decrease of the ex- 
citon occupation for high pulse areas. 

For all chirp values, we find that for a biexciton shift 
of Eb = —4 meV the four-level model agrees with the 
two-level model. This is because a biexciton shift of 
Eb = —4 meV leads to an energy separation above the 
cut-off frequency of the carrier-phonon coupling and thus 
no phonon-assisted biexciton generation is possible. The 
same effect takes place for the positive biexciton shift. The 
results for Eb = 2 meV (not shown in the pictures) show 
perfectly agreement with the case of circular polarization. 

4 Summary In summary, we have investigated the 
exciton dynamics in a QD adiabatically controlled by 
a chirped laser pulse. For circular polarization we have 
compared results obtained from three different simulation 
methods namely (i) TCL, (ii) CE and (iii) PI. The results 
obtained from CE and PI show very good agreement, even 
at T = 50 K. The TCL agrees perfectly with the other 
two methods at low temperatures and agrees reasonably 
at moderate temperatures, but at high temperatures a mis- 
match of up to 30% was found. The agreement between 
these methods is better for the ARP dynamics than for the 
Rabi oscillations. We have also studied exciton dynamics 
for a linearly polarized laser pulse. We have shown, that 
new decoherence paths related to phonon-assisted biexci- 
ton generation open up, which are efficient for biexciton 
shifts of the order of —2 meV. For higher biexciton shifts 
the influence of the biexciton becomes negligible and the 
system can be reduced to a two-level model. 
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tion and by the TEAM programme of the Foundation for 
Polish Science, cofinanced from the European Regional 
Development Fund. 
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